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CHAPTER  1 


INTRODUCTION 


Design  of  a  system  to  meet  a  demand  or  need  is  a  function  which  must 
be  addressed  in  its  entirety.   The  total  system  must  be  examined  and  the 
interplay  of  the  design  variables  and  system  parameters  must  be  intimate- 
ly considered.   This  systems  approach  is  not  the  usual  case  today.   The 
system  approach  is  often  attempted  after  the  system  is  ready  for  delivery 
or  for  redesign. 

In  major  systems  the  piece-meal  design  procedures  used  are  often 
the  cause  of  inadequate  system  performance.   The  systems  and  units  often 
cannot  be  economically  maintained  or  operated  in  the  real  world  environ- 
ment.  The  logistics  and  support  for  the  system  is  poor  and  generally 
determined  after  the  system  comes  into  being. 

Performance,  cost,  and  supportability  can  be  considered  key  ingre- 
dients in  the  development  of  any  system.   Optimization  of  these  factors 
simultaneously  is  a  difficult  task.   Trade-offs  are  always  required. 
There  exists  a  need  for  more  interplay  between  design,  system  analysis, 
and  supportability.   The  goals  established  in  the  areas  of  performance 
requirements  in  terms  of  effectiveness,  cost,  and  supportability  must 
initially  be  addressed  in  the  design  process.  Design  and  maintainability 
engineering  must  translate  goals  into  effective  means  to  make  the  con- 
cepts viable.   This  is  an  interative  time  consuming  process,  but  impor- 
tant as  a  first  step. 


This  study  developed  a  model  and  algorithm  to  be  utilized  as  a 
basis  for  optimal  design  of  a  finite  population  queuing  system  deployed 
to  meet  a  demand.   The  sensitivity  and  constraints  of  the  system  para- 
meters and  design  variables  are  then  investigated  with  the  algorithm. 

Operational  System 

The  system  under  study  is  illustrated  in  Figure  1,  and  is  designed 
to  meet  the  demand,  D.  Units  in  the  system  are  in  either  of  two  cate- 
gories; they  are  available  operationally  (shaded  area)  or  they  are  in  a 
repair  state  (in  queue  or  service) .   In  this  model  the  units  are  all 
repairable  when  they  fail.   It  is  assumed  that  there  are  no  units  beyond 
economical  repair  or  units  destroyed. 

All  units  which  fail  enter  either  the  intermediate  repair  level 
(IRL)  or  the  depot  repair  level  (DRL)  and  will  be  returned  to  operational 
use.   This  model  will  determine  the  optimal  number  of  repair  channels  at 
both  the  IRL  and  DRL,  the  optimum  number  of  units  which  must  be  provided 
to  meet  the  demand  utilizing  the  two  level  maintenance  concept,  and  the 
optimum  replacement  age.   If  it  is  calculated  to  be  more  economical  to 
have  no  depot  level  but  only  an  IRL,  the  output  of  the  algorithm  will 
specify  this  alternative.   The  model  is  flexible  in  that  allowances  are 
made  to  permit  a  change  of  the  pertinent  parameters  in  order  to  determine 
the  optimum  design  for  specific  situations. 

One  of  the  model's  primary  outputs  is  the  determination  of  the 
optimum  throw-away  age.   In  Figure  1,  new  units  enter  the  system  and 
remain  until  the  nth  year.   The  number  of  units  is  assumed  to  be  a 
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constant  and  the  number  in  each  age  is  also  a  constant  equal  to  the 
ratio  of  the  total  population  size  and  the  optimum  throw-away  age. 

A  system  approach  must  be  utilized  to  analyze  and  integrate  all 
the  various  components  and  variables  involved.   In  the  sensitivity  anal- 
ysis of  the  model,  the  effects  of  mean-tirae-between-failure  (MTBF) , 
mean-time-to-repair  (MTTR),  percent  allocation  for  the  IRL  and  DRL, 
depreciation  mode,  interest  rate,  and  depot  MTTR  will  be  investigated. 
In  the  constraint  analysis,  shorter  desired  life  for  the  units  and 
population  size  variations  will  be  analyzed. 

The  system  approximates  the  operation  of  many  systems  in  existance. 
An  automobile  rental  service,  numbers  of  aircraft  in  a  fleet  necessary 
to  meet  operational  missions,  or  industrial  equipment  required  to  pro- 
duce specified  outputs  are  a  few  examples.   This  model  is  flexible 
enough  to  fit  many  systems  with  the  basic  overall  characteristics  de- 
scribed above. 

Scope  and  Assumptions 

The  model  deals  with  the  design  variables  and  system  parameters 
which  have  an  impact  on  the  ability  of  the  system  to  meet  the  demand. 
However,  the  system  transient  phase  is  not  considered,  this  is  a  steady 
state  model.   The  portion  of  the  life-cycle  to  be  considered  is  the 
steady  state  operation  which  occurs  between  the  introduction  phase  and 
the  phase-out  phase  of  the  system. 

Other  assumptions  which  are  made  concerning  the  system  are  as 
follows: 


1)  The  number  of  units  deployed  to  meet  the  demand  is 
not  constrained. 

2)  The  population  may  be  made  up  of  units  of  many  ages 
but  are  otherwise  homogenous. 

3)  All  failed  units  are  repairable,  no  failed  units 
will  be  discarded  from  the  system. 

A)  MTBF  and  MTTR  change  constantly  for  a  unit,  however 
they  are  assumed  to  be  constant  over  each  yearly 
period,  the  figures  used  are  yearly  averages. 

5)  The  population  is  considered  finite  for  calculation 
purposes,  and  will  be  held  below  50  units  due  to 
computer  limitations. 

6)  Spare  units  will  be  in  the  maintenance  system  for 
training  but  will  not  be  considered  as  a  part  of  the 
total  population  available  to  meet  the  demand. 

System  Effectiveness 

The  mathematical  model  formulated  in  this  study  has  as  its  basic 
goal  the  determination  of  a  system  configuration  which  will  effectively 
meet  the  demand.   The  measure  of  its  ability  to  meet  this  demand 
effectively  is  the  cost  incurred.   The  cost  equations  will  be  based  on 
the  annual  equivalent  cost  concepts  of  engineering  economics. 

There  are  two  classes  of  variables  involved  in  this  system,  those 
under  direct  control  of  the  designer  and  those  that  are  not  directly 
under  his  control.   Both  types  are  of  primary  concern  in  this  development 


Those  which  are  directly  under  the  designers  control  are  design  variables 
whereas  those  not  under  his  direct  control  are  system  parameters.   System 
effectiveness  depends  directly  on  the  values  chosen  for  the  design 
variables  in  the  face  of  system  parameters. 

System  Parameters 

The  eleven  variables  not  directly  under  the  control  of  the  designer 
are  (1)  initial  cost  of  the  units,  (2)  demand,  (3)  initial  facility 
cost,  (4)  MTBF,  (5)  MTTR,  (6)  percent  allocation  of  failed  units  to  IRL 
and  DRL,  (7)  depreciation  method,  (8)  depot  MTTR,  (9)  interest  rate, 
(10)  service  distribution,  and  (11)  the  failure  distribution. 

The  cost  of  spare  parts,  training,  documentation,  spare  units  for 
training  are  also  system  parameters  in  this  model.   The  cost  of  unmet 
demands  is  a  function  of  the  interaction  of  these  parameters  and  the 
quality  and  quantity  of  maintenance  and  spares. 

The  MTBF  and  MTTR  are  not  directly  controlled  by  the  designer, 
however  he  may  decide  to  have  them  varied  by  paying  the  appropriate 
costs  of  redesign  and  engineering  to  either  improve  or  degrade  avail- 
ability.  This  study  does  not  take  into  account  scheduled  maintenance. 
The  MTTR  will  include  the  transit  time  and  repair  time  but  not  the 
waiting  time  in  the  queue. 

Design  Variables 

The  four  design  variables  are  (1)  throw-away  age  of  the  units,  (2) 
population  size,  (3)  the  number  of  IRL  channels,  and  (4)  the  number,  of 


DRL  channels.   The  number  of  units  within  an  age  group  will  be  assumed 
equal.   For  example,  a  system  containing  a  population  of  units  with  a 
maximum  age  of  three  years  comes  into  existence  by  procuring  one  third 
of  the  units  in  each  of  three  consecutive  years  and  replacing  one  third 
each  year  thereafter.  Under  this  assumption  a  population  exists  which 
will  exert  a  steady  load  on  the  maintenance  facilities.   It  will  be 
assumed  that  the  system  has  reached  its  steady  state  by  some  prior  pro- 
curement action. 

The  number  of  service  channels  to  be  determined  for  both  IRL  and 
DRL  is  of  critical  importance.   The  severity  of  the  failure  will  determine 
which  level  of  maintenance  the  unit  will  enter;  the  channels  at  the  IRL 
and  DRL  will  have  their  own  repair  times.   Each  channel  will  have  equal 
capability  so  that  the  repair  time  is  equivalent  regardless  of  the 
type  of  failure.   The  arrivals  will  be  as  determined  by  the  proportion 
of  failures  relegated  to  each  level.   The  arrivals  will  enter  queues  at 
the  appropriate  facility  and  be  accepted  on  a  first  come  first  served 
basis.   The  channels  required  in  either  IRL  or  DRL  will  be  parallel 
service  channels  and  they  can  serve  units  simultaneously.  Demand  is  not 
a  design  variable  as  it  is  the  driving  force  for  the  development  of  the 
system.   It  is  the  requirement  for  which  the  system  is  being  designed 
to  meet. 

Survey  of  the  Literature 

The  most  significant  documents  dealing  with  the  procurement  and 
acquisition  of  units  or  entire  systems  have  been  written  for  the   , 


Department  of  Defense.   They  are  "Life  Cycle  Costing  Procurement  Guide" 
(4),  and  "Life  Cycle  Costing  Guide  for  System  Acquisition"  (5).   The 
former  guide  is  designed  to  aid  in  the  procurement  of  the  units  and 
materials  which  make  up  only  components  of  the  total  system;  the  latter 
is  designed  to  apply  the  life  cycle  cost  concept  to  the  entire  system. 
These  models  are  applicable  to  decision  making  analysis;  i.e.,  comparative 
cost  analysis  of  alternative  systems,  and  analysis  of  the  same  system 
under  various  support  or  replenishment  concepts.   These  models 
attempt  to  determine  the  optimum  replacement  age  for  the  systems  not 
the  units  in  the  system.   The  number  of  units  to  be  acquired  must  be 
specified  in  both  of  these  models. 

There  are  a  number  of  simulation  models  available  to  investigate 
life  cycle  costs  of  a  system.   These  models  attempt  to  find  the  optimum 
number  of  repair  facilities,  number  of  spares  required,  and  personnel 
required  to  support  a  specified  system.   The  "System  Cost  and  Operational 
Resources  Evaluation  (SCORE)"  model  developed  by  the  Naval  Air  Develop- 
ment Center  (12),  is  an  example  of  such  a  model.   It  estimates  the  life 
cycle  costs  of  research  and  development,  test  and  evaluation,  investment 
and  operations  for  fifteen  years  for  the  various  component  elements, 
and  aggregates  these  into  a  total  cost  estimate  for  the  system.   The 
models  mentioned  above  require  large  amounts  of  data  to  adequately  per- 
form their  analysis  and  require  extremely  long  computer  runs. 

Much  of  the  effort  in  the  area  of  system  costs  has  been  directed 
to  defining  the  reliability  and  maintainability  aspects  of  the  systems. 
Armstadter  (1) ,  and  Gnedenko  (8)  do  mathematical  analysis  of  systems 


to  attain  the  probable  best  combinations  of  reliability  and  maintain- 
ability but  do  not  consider  any  cost  functions  and  their  interrelation- 
ships to  these  parameters  and  effects  on  the  system  performance.   They 
define  MTBF  as  a  measure  of  reliability  and  MTTR  as  a  measure  of  main- 
tainability.  Seller  (15)  does  an  analysis  of  a  system  assuming  that  the 
probability  functions  of  cost  are  triangular  in  shape  and  that  there  is 
only  one  optimal  cost  with  all  other  costs  linearly  distributed  on 
either  side  of  the  optimum.   Briggs  (3)  did  work  in  determining  the 
reliability  of  systems  and  its  effect  of  cost.  His  effort  is  in  areas 
which  are  generally  post  design  phases,  after  the  system  has  been  in 
operation  for  a  period  of  time. 

MILHANDBOOK  472  (11)  defines  the  various  distributions  of  failure 
and  repair  times  and  rates.   Specific  statistical  distributions  are 
listed  which  correspond  to  those  used  in  the  model.   The  basis  for  the 
queuing  equations  used  in  this  model  is  similar  to  that  described  and 
utilized  by  Jelen  (9)  and  also  those  problems  in  the  area  of  queuing 
as  discussed  by  Kendall  (10) . 

Wortham  (17)  bases  his  analysis  of  system  cost  on  the  initial  cost 
of  the  system  and  support  costs.  In  his  analysis,  the  availability  of 
the  individual  units  was  considered  as  not  critical  when  the  value  was 
between  80  and  85  percent  because  then  there  would  exist  ample  time  to 
perform  maintenance  and  support  functions.  However  if  the  availability 
was  necessarily  between  95  and  97  percent,  then  the  design  reliability 
and  maintainability  were  considered  very  critical. 
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In  Goldman  (7),  the  use  of  a  particular  definition  as  to  what  are 
support  costs  can  change  the  various  component  life  cycle  cost  calcula- 
tions considerably.   Many  large  systems  support  costs  often  exceed  the 
cost  of  the  equipment  by  more  than  ten  times,  even  100  times  is  not  rare, 
Also,  Brandt  (2)  developed  a  model  which  is  utilized  to  predict  the 
down-time  of  a  unit  as  a  function  of  its  utilization,  and  whether  it  is 
utilized  fully  or  partially  when  operated. 

In  summary,  a  review  of  the  available  literature  on  subjects 
related  to  this  thesis  does  not  reveal  any  method  of  optimizing  all 
the  design  variables  for  a  general  operating  system  as  addressed  herein. 
Further,  none  of  the  analyses  consider  the  level  of  maintenance  concept 
or  the  impact  of  the  various  variables  on  system  effectiveness. 


CHAPTER  2 
MATHEMATICAL  MODEL  AND  COMPUTER  ALGORITHM 

In  attaining  the  primary  goal  of  this  study  it  was  necessary  to 
develop  a  model  and  an  associated  algorithm  which  would  permit  determi- 
nation of  the  throw-away  age,  population  size,  and  number  of  service 
channels  in  IRL  and  DRL  to  meet  demand  in  an  optimum  manner.   The  model 
is  cost  based;  it  incorporates  the  annual  equivalent  cost  (AEC)  of  the 
maintenance  facilities  and  the  annual  equivalent  cost  of  the  units  plus 
the  annual  cost  for  non-fulfillment  of  demands. 

The  following  assumptions  are  made  in  the  development  of  the 
mathematical  model: 

1)  Inter-arrival  times  are  exponentially  distributed  at 
both  repair  levels. 

2)  The  repair  times  are  exponentially  distributed  at 
both  repair  levels. 

3)  The  inter-arrival  times  are  statistically  independent 
of  the  arrival  times,  at  both  repair  levels. 

4)  The  service  channels,  at  both  IRL  and  DRL  are  parallel 
and  each  is  capable  of  performing  the  maintenance 
activities  associated  with  its  respective  area  of 
responsibility. 

5)  The  population  size  will  always  be  larger  than  the 
number  of  service  channels. 
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6)  Each  channel  performs  service  on  one  unit  at  a  time. 

7)  Units  are  serviced  in  the  repair  channels  in  a 
first  come  first  served  mode. 

8)  Units  completing  repair  return  to  operation  with  the 
same  operational,  reliability  and  maintainability 
characteristics  established  for  its  associated  age. 

9)  The  cost  of  the  repair  facility  varies  linearly  with 
the  number  of  channels.   The  initial  investment  cost 
in  land  and  associated  improvements  is  constant  for 
up  to  three  channels;  the  base  cost  of  the  facilities 
is  increased  for  each  additional  increment  of  three 
channels,  this  pertains  to  both  levels  of  maintenance. 

10)  The  out-of -operation  cost  is  a  linear  function  of  the 
number  of  units  not  available  to  meet  demand . 

11)  The  administrative,  labor,  and  operating  expenses  for 
the  IRL  and  DRL  facilities  increase  linearly  with  the 
number  of  channels  being  required. 

12)  The  cost  of  training,  initial  spares,  and  recurring 
support  are  included  as  proportions  of  the  total  annual 
equivalent  cost  (TUAEC)  and  the  IRL  facility  costs. 

Population  Annual  Equivalent  Cost 

To  determine  the  cost  of  ownership  of  the  population  of  units  on 
an  annual  equivalent  cost  basis,  the  fundamental  equations  of  engi- 
neering economy  are  used  with  the  notation  modified  for  this  thesis  (16). 
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TUAEC  ■  Population  annual  equivalent  cost 

P  =  First  or  acquisition  cost  of  the  unit 

SV  =  Estimated  salvage  value  of  the  unit 

BT  =  Book  value  at  end  of  year  AGE 

AGE  =  Year  of  oldest  throw-away  age,  AGE  <_  NP 

NP  =  Estimated  life  of  the  unit 

N  =  Number  of  units  in  the  population 

i  ■  Interest  rate 


TUAEC  =  <(P-BT) 


i(H-i)AGE 
(l+i)AGE  -1 


+  BT(i)  >  N 


The  book  value  of  the  units  at  the  time  of  their  disposal  will 
depend  upon  the  above  equation  and  the  method  of  depreciation  used. 
The  methods  of  calculating  the  depreciation  of  an  item  over  time  are 
given  in  the  following  sections. 

Straight  line  method.   The  straight  line  method  of  depreciation 
assumes  that  the  value  of  an  asset  decreases  at  a  constant  rate. 


BT  =  P  -  AGE  £=§¥■) 


Sum  of  digits  method.   The  sum  of  the  years  digits  depreciation 
method  assumes  that  the  value  of  an  asset  decreases  at  a  decreasing 
rate. 


bt  -  (P-sv)  (^)  (Sgp.)  ♦  SV 


14 


Sinking  fund  method.   The  sinking  fund  method  of  depreciation 
assumes  that  the  value  of  an  asset  decreases  at  an  increasing  rate, 


BT  =  P  -  (P-SV) 


(l+i)^  -1 


(l+i)AGE  -1 


Zero  book  value.  Although  not  specifically  mentioned  in  the 
references,  in  some  circumstances  it  is  desirable  to  set  a  zero  value 
on  an  asset.   This  will  take  into  account  the  situation  where  units 
are  retired  (at  any  age)  with  no  appreciable  value.   From  the  above 
equations,  the  total  annual  equivalent  cost  of  all  units  can  be 
calculated  knowing  the  age  of  the  oldest  units. 


Repair  Facility  Annual  Equivalent  Cost 

The  annual  equivalent  cost  for  both  IRL  and  DRL  facilities  must 
be  considered.  The  basic  formula  and  the  notation  to  be  used  are  as 
follows: 

LFS 


FAEC  =  (PFAC-FFAC) 


Kl+1) 


(1+i) 


LFS 


+  FFAC(i) 


where, 

FAEC 

PFAC 

FFAC 

LFS 

i 


facility  annual  equivalent  cost  per  channel 

first  cost  of  the  facility  per  channel 

estimated  salvage  value  of  the  facility  per  channel 

estimated  life  of  the  facility 

interest  rate 
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This  equation  is  similar  to  the  one  derived  previously  for  TUAEC.  This 
cost  is  a  linear  function  of  the  number  of  repair  channels  and  has  a 
non-zero,  positive  value  when  the  number  of  repair  channels  are  cal- 
culated.  The  value  of  the  depot  facility  annual  equivalent  cost  (FAECD) 
for  the  situation  where  there  are  zero  depot  channels  is  zero,  there 
will  be  no  investment  of  any  type  at  that  level.  However,  the  cost  of 
the  IRL  will  be  increased  in  the  algorithm  to  take  into  account  the 
increased  repair  requirements  at  that  level  as  its  maintenance  respon- 
sibility would  be  much  greater. 

In  the  calculation  of  the  facility  cost  the  various  types  of 
equipment  and  buildings  must  be  taken  into  consideration,  these 
different  values  and  different  life  expectancies  are  handled  as  a 
summation  for  each  level  of  repair.   Each  facility  can  handle  only  a 
given  number  of  channels.  A  factor  is  added  to  the  algorithm  which 
will  increase  cost  as  the  number  of  channels  is  increased.   This  cost 
is  a  function  of  the  number  of  channels  for  both  types  of  facilities. 

Let 

FAECN  =  total  facility  annual  equivalent  cost 

DY    =  yearly  administrative,  manpower,  and  overhead 

costs  per  channel 
CC    =  number  of  IRL  channels;   CC  >  0 

(for  DRL,  DC  =  number  of  DRL  channels;  0  £  DC  £  6) 
ALEO  =  first  cost  of  establishing  facilities 

Then, 

FAECN  =  CC(FAEC  +  DY)   +  ALEO 
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This  equation  is  utilized  with  two  sets  of  notation  in  the  model 
to  facilitate  handling  of  IRL  and  DRL.   The  notation  for  FAEC  which 
contains  a  D,  pertains  to  the  depot  facilities. 

Out-of-Operation  Cost 

When  units  are  not  available  to  meet  a  specified  demand  some  cost 
is  incurred.   The  number  of  units  available  at  any  given  time  is  a 
random  variable,  but  is  of  unknown  distribution.   The  out-of -operation 
cost  is  therefore  a  function  of  the  demand  and  the  number  of  units 
available.   This  number  is  a  function  of  the  failure  rate,  the  time  to 
repair,  the  population  size,  the  number  of  unit  age  groups,  the  expect- 
ed ratio  of  major  and  minor  failures,  and  the  number  of  IRL  and  DRL 
facilities  being  provided. 

The  first  parameters  which  need  to  be  determined  are  the  MTBF  and 
MTTR  values  at  both  levels  of  maintenance,  for  the  populations  at  each 
level.   It  is  assumed  that  the  value  of  these  parameters  is  known  or 
can  be  estimated  initially.   They  are  input  system  parameters  and 
extend  over  the  life  of  the  units.   These  inputs  are  on  a  yearly  basis. 

Since  it  has  been  assumed  that  the  number  of  failures  per  unit 

time  and  the  repairs  per  unit  time  are  Poisson  distributed,  a  theorem 

by  Parzen  (13)  can  be  used  for  the  cumulative  sums  for  the  various 

unit  ages.   The  theorem  states  that  if  S  =  Xn+X0+...+X  is  the  sum  of 
to  n    1  2      n 

n  independent  random  variables,  if  k  =  l,...,n;  X,  is  Poisson  distri- 
buted with  parameter  A.+...+X  .  Also  the  mean  of  S  is  equal  to  the 
sum  of  the  individual  means. 
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Let 

MTBF.  =  mean  time  between  failures  for  the  jth  unit  age 

MTTR  ■  mean  time  to  repair  value  for  the  jth  unit  age 

n  =  number  in  the  unit  age 

N  =  population  size 

Assuming  each  throw-away  age  contains  the  same  number  of  units, 
the  calculations  necessary  to  find  the  MTBF  and  MTTR  are  as  follows: 


n    MTBF.  (-) 
MTBF  =  I         L-S-  (2.1) 


n   MTTR.  (-) 
MTTR.  =   ^   J— -  (2.2) 


Once  these  values  are  determined,  the  expected  number  of  units 
operational  can  be  determined  by  using  the  Finite  Queuing  Tables  of 
Peck  and  Hazlewood  (14)  or  computed  using  the  equations  developed  in 
Appendix  A.   The  results  derived  in  Appendix  A  are  used  for  this 
algorithm  and  the  notation  given  below  appears  in  the  model. 


f     =  MTBF 


-  =  MTTR 


x 

p  =  - 

V 


W  =  expected  time  a  unit  waits  in  a  queue 


18 


W  = 


N-l       ,- 

I       pi££±l)  N_n 

L        K     cc  * 
n=cc 


N! 


ccCC(^-) 


_(N-n)  !J    cc ! 


cc-1 


N! 


T    (N-n)   .J"1  ,,(p)n  +  I        (N-n)  ,._M \, 
_0       n!(N-n)I  K     «_««       (N-n)! 


cc 


cc 


© 


n 


n=cc 


ccl 


(2.3) 


L  =  expected  number  of  units  in  a  repair  level  (queue 
plus  repair  channel) 


L  =  N  (  1  - 


(1/y+l/A) 
(1/u+lA+W  ) 


1  - 


l/v 


(1/y+l/X) 


(2.4) 


The  expected  number  of  units  available  is  determined  using  the 
following  equation: 

TJ  =  XN  -  (L  +  LL) 


where , 
TJ 
XN 
L 
LL 


expected  number  of  units  available 
total  unit  population 
expected  number  of  units  in  IRL 
expected  number  of  units  in  DRL 


This  number  is  utilized  in  calculating  the  cost  incurred  for  not  meeting 
a  demand.   Since  the  number  of  units  available  is  a  random  variable,  it 
has  some  type  of  distribution  (which  is  not  determined  herein) .   It 
is  assumed  that  there  is  a  degree  of  variability  in  the  number  of  units 
available  which  could  alter  the  cost  of  not  meeting  the  demand.  For 
simplicity  this  cost  will  be  assumed  to  be  a  linear  function  of  th'e 
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number  of  units  demanded  and  not  available. 


Let 


C   ■  cost  of  unmet  unit  demands  per  day 
Z   =  cost  of  unmet  demand  per  unit  time  period 
D   =  expected  demand 

Q   =  degree  of  variability  of  the  number  of  units 
available  at  any  given  time 


C   =  < 

w 


(Q  •  Zn)  -  Z  (TJ  -  D) ,   Q  >  (TJ  -  D) 


,   Q  <  (TJ  -  D) 


The  system  cost  includes  the  training  of  personnel  to  operate  and 
maintain  the  units  which  are  procured  for  the  system.   Seven  percent 
of  the  total  unit  and  facility  annual  equivalent  cost  is  added  to  the 
system  cost  for  training  and  initial  support.  The  spare  parts  and 
piece  parts  are  a  recurring  expense,  twenty  percent  of  the  TAEC  of  the 
units  and  facilities  is  used  in  this  model  as  their  cost  estimate. 
The  grand  total  annual  equivalent  cost  is  the  sum  of  all  costs  pre- 
viously specified  and  is  expressed  as  follows: 

TTAECN  =  TTAEST  +  FFAEST  +  STRN  +  RSUPN  +  WCN 

where 

TTAECN  =  grand  total  annual  equivalent  cost 

TTAEST  =  total  annual  equivalent  cost  of  intermediate 

r 

plus  depot  level  units 


20 


FFAEST  =  total  annual  equivalent  cost  of  the  intermediate 

plus  depot  facilities 
STRN   =  training  and  initial  support  cost 
RSUPN  =  recurring  support  cost 
WCN    =  cost  due  to  not  meeting  a  demand 

The  Computer  Algorithm 

The  program  coding  which  was  developed  for  this  algorithm  is 
exhibited  in  Appendix  B.l.   The  program  is  made  up  of  a  main  routine 
and  three  subroutines.   Each  of  these  areas  will  be  described  in  the 
following  sections. 

Main  program.   The  primary  objective  of  this  program  is  to  minimize 
(optimize)  the  TAEC.   The  goal  of  the  main  program  is  to  execute  a 
search  technique  which  will  optimize  the  total  cost  function  over  the 
range  of  variables  delineated  by  the  designer,  and  to  act  as  the  handler 
of  the  other  routines.   The  total  system  cost  function  is  unimodal. 
There  is  only  one  true  minimum  for  any  given  set  of  system  parameters. 
The  number  of  variables  over  which  the  search  must  progress  is  four: 
the  number  of  IRL  channels,  the  number  of  DRL  channels,  total  popu- 
lation size,  and  the  throw-away  age  of  the  units. 

To  commence  the  search  it  is  necessary  to  determine  the  service 
factors  to  be  used  in  the  IRL  and  DRL.  This  equation  was  originally 
developed  and  used  by  Peck  and  Hazlewood  (14) . 
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Let 


X  =  Service  factor  (IRL)   = 


XMTTR  +  XMTBF 


XD  =  Service  factor  (DRL)   = 


DMTTR  +  DMTBF 


Using  the  number  of  channels  in  IRL  as  the  starting  function  for 
the  search  routine,  the  number  of  channels  from  which  to  commence  the 
search  is  calculated  as  follows: 

CC  =  (N  •  X)  +  4  [N  •  X(l-X)]1/2 

The  search  is  run  through  the  range  of  IRL  channels  looking  for 
the  minimum  TAEC.  During  this  search  the  number  of  DRL  channels  is 
maintained  at  a  constant  value.  After  each  minimum  TAEC  is  determined 
for  a  population  size,  the  DRL  channels  available  is  incremented.  A 
run  is  made  with  each  succeeding  DRL  value.   Only  the  minimum  set  of 
TAEC  data  for  each  combination  of  IRL  and  DRL  channels  is  stored  in 
the  computer.   This  minimum  can  be  visualized  as  a  point  on  a  plane  in 
three  dimensional  space,  the  plane  being  for  a  given  population  size 
and  throw-away  age.   For  each  throw-away  age  considered  there  will  be 
a  set  of  outputs;  these  must  be  visually  scanned  for  the  minimum  TAEC 
for  that  throw-away  age.  When  all  population  sizes  and  throw-away  ages 
which  are  being  considered  have  been  searched,  the  minima  must  be 
checked  again  to  identify  the  fourth  minimum;  this  is  the  ultimate 
objective  of  the  main  program,  the  global  minimum  TAEC.   The  results 
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indicate  the  number  of  IRL  AND  DRL  channels,  and  the  associated  cost 
of  that  population  size  and  throw-away  age. 

Subroutine  RAIDS.  This  subroutine  determines  the  rate  of  arrivals 
and  service  at  IRL  and  DRL.   It  provides  the  aggregate  XMTBF  and 
XMTTR,  DMTBF  and  DMTTR  values  for  IRL  and  DRL  for  each  population  size 
and  throw-away  age  combination  each  time  they  are  required  in  the  com- 
putations.  The  equations  in  (2.1)  and  (2.2)  are  utilized  in  the 
execution  of  this  subroutine.   This  routine  includes  the  calculation 
of  the  inherent  availability  of  the  units  in  the  intermediate  level. 
This  is  the  designed  availability  (not  the  actual  availability)  which 
is  determined  after  the  system  has  been  in  use  for  some  period  of  time. 
The  availability  is  contained  in  the  printout  for  each  population  size. 
It  is  computed  as  follows: 

XAI  =     3MrBF 


XMTBF  +  XMTTR 


XAI  =  Inherent  availability 

Subroutine  UNAVAIL.   This  subroutine  determines  the  number  of 
units  available  from  a  given  population  size  utilizing  the  appropriate 
values  determined  in  RAIDS.   The  important  result  of  the  calculations 
using  equations  (2.3)  and  (2.4)  is  the  expected  number  of  units  not 
available.   This  subroutine  is  limited  in  the  population  size  it  can 
accommodate  due  to  the  size  of  the  factorials  which  must  be  calculated; 
55  factorial  is  the  limit  for  the  computer  used  with  this  model.  A 
limit  may  also  be  reached  if  the  desired  availability  of  the  units  is 
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high.   Generally  the  program  will  function  adequately  when  the  ratio 
of  MTTR  to  MTTR  plus  MTBF  is  greater  than  two  tenths.   The  maximum 
population  size  which  can  be  run  on  the  IBM-370  Computer  System  irre- 
spective of  the  ratio  described  above  is  43  units. 

Subroutine  INVEST,  All  cost  factors  mentioned  previously  are 
included  in  the  overall  system  cost  and  are  calculated  and  summed  in 
this  routine.   It  includes  the  MOD  designation  for  the  book  value 
calculations.  This  routine  calculates  all  the  costs  for  the  IRL  and 
DRL  facilities  which  includes  an  escalation  factor  for  the  increased 
number  of  channels  in  either  facility.  The  TUAEC  for  the  units  in  both 
levels  of  maintenance,  training  and  initial  support,  and  recurring 
support  costs  are  also  computed  in  this  routine.  All  costs  are  totaled 
and  the  result  is  the  TAEC  for  the  particular  combination  of  para- 
meters being  investigated.   The  basic  time  unit  used  throughout  all 
the  routines  is  one  day  in  a  365  day  year. 


CHAPTER  3 

STANDARD  SYSTEM  DESIGN  EXAMPLE 

This  chapter  presents  a  standard  system  design  formulated  to 
utilize  the  mathematical  model  and  exercise  the  optimization  algorithm 
developed  in  the  previous  chapter.   The  standard  system  parameter  in- 
puts must  be  obtained  by  the  designer  and  represent  a  system  with  which 
he  is  concerned.  He  will  receive  from  the  programmed  algorithm  the 
optimal  combination  of  IRL  and  DRL  channels,  the  minimum  TAEC,  and  when 
to  throw-away  the  units  for  the  most  effective  and  efficient  system 
utilization. 

System  Parameter  Input  Values 

Appendix  B.2  gives  the  standard  system  design  input  values.  The 
demand  is  20  units  per  day.   The  maximum  age  of  the  units  in  the  system 
is  15  years.   The  depreciation  code  is  two  (straight  line  method).   The 
number  of  units  to  be  repaired  at  each  level  is  allocated  75  percent  to 
IRL  and  25  percent  to  DRL. 

The  depot  MTTR  is  2.5  days  and  is  constant  over  the  life  of  the 
units.  All  units,  when  sent  to  depot,  will  be  overhauled  and  returned 
to  service  with  the  same  performance  and  maintainability  character- 
istics as  their  original  age  group.   The  depot  MTBF  will  be  the  same 
as  the  IRL  MTBF;  all  units  will  have  the  same  reliability. 

The  severity  of  the  failures  is  not  a  function  of  the  MTBF  but 
is  a  random  factor  in  the  unit  design  which  is  unpredictable.  The' 
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value  of  MTBF  and  MTTR  for  the  IRL  are  given  next.  A  value  for  each 
year  is  estimated  and  is  the  average  value  for  the  units  in  that  throw- 
away  age,  i.e.,  in  the  second  year  the  MTBF  is  6.5  days  and  the  MTTR 
is  1.3  days. 

The  MTBF  of  the  units  in  this  system  are  assumed  to  have  the  "bath 
tub"  failure  hazard  characteristic.  During  the  first  few  years  the 
MTBF  will  improve  and  reach  its  designed  value.   It  will  slowly  degrade 
as  the  units  age.  As  the  units  near  the  end  of  their  projected  life 
the  failure  rate  will  increase  rapidly. 

The  MTTR  will  also  vary.  At  initial  introduction  the  time  to 
repair  will  be  longer  than  desired  for  the  units  and  will  reach  its 
design  value  about  the  same  time  that  the  MTBF  achieves  its  design 
value.   The  MTTR  will  increase  as  the  units  age  due  to  in-depth  repairs 
becoming  more  prevalent. 

The  cost  data  used  in  this  example  will  be  as  follows:   The 
initial  unit  cost  is  $25,000,  with  a  maximum  life  of  15  years.  The 
facility  costs  are  separated  into  three  categories  for  both  the  IRL 
and  DRL.   The  first  of  each  set  is  for  buildings  and  related  equipment, 
with  the  IRL  expected  to  last  40  years  and  the  DRL  30  years.   The 
large  repair  equipment,  test  benches,  and  test  equipment  is  the  second 
set  and  will  be  usable  for  10  and  8  years  respectively.   The  third  set 
is  for  tools,  test  kits,  and  small  support  equipment  and  will  last  5 
years  at  both  levels.  All  facility  costs  are  listed  in  Appendix  B.2 
and  are  calculated  in  the  subroutine  INVEST  on  a  per  channel  basis. 
The  interest  rate  on  all  investments  is  10  percent. 
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The  expected  yearly  administrative  labor  and  overhead  costs  are 
$25,000  and  $62,500  respectively  at  IRL  and  DRL.  The  cost  of  land, 
utilities,  etc.  are  included  as  inputs  and  are  computed  in  the  algorithm 
as  a  step  function  dependent  on  multiples  of  three  repair  channels. 
The  initially  assumed  cost  is  $25,000  and  $75,000  respectively. 

The  cost  for  not  meeting  a  demand  is  $200  per  unit  per  day.   The 
units  available  will  vary  by  a  factor  of  four  when  there  is  no  DRL 
and  by  a  factor  of  2.5  when  both  levels  are  available.   If  the  number 
of  available  units  is  equal  to  or  greater  than  22.5  there  will  be  no 
cost  incurred,  however  should  the  number  of  units  available  be  less 
than  22.5,  a  cost  will  be  incurred. 

Design  Variable  Output  Results 

Appendix  B.3  gives  the  results  using  the  algorithm  and  standard 
system  design  example  input.   The  printout  shows  all  the  costs  which 
make  up  the  minimum  cost  combination  of  IRL  and  DRL  channels  for  a 
given  throw-away  age  and  population  size.  To  allow  for  better  under- 
standing of  the  output  of  the  model,  the  seventh  throw-away  age 
listing  and  a  population  of  34  units  will  be  discussed  in  detail,  then 
the  entire  throw-away  age  range  is  analyzed. 

The  algorithm  starts  at  throw-away  age  one  and  proceeds  to  the 
maximum  age  designated  in  the  input  of  15  years.   The  unit  inherent 
availability  for  this  case  is  86.21  percent.   The  number  of  IRL 
channels  is  5  and  the  DRL  channels  is  1.   The  queuing  routine  described 
in  Chapter  2  is  used  to  determine  the  number  of  units  that  will  no,t 
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be  available.   The  channels  are  searched  until  a  minimum  TAEC  is  found. 
For  this  problem  the  cost  computed  for  IRL  is  as  follows:   TUAEC  is 
$106,083.81,  and  the  FAEC  is  $176,957.88.   For  the  DRL,  the  TUAEC  is 
$35,361.28  and  the  FAEC  is  $15,006.44.   The  training  and  initial  support 
cost  is  $19,812.91,  while  the  recurring  support  cost  is  $56,608.33.   The 
cost  for  not  meeting  the  demand  is  $0.00.  Also  listed  are  the  number 
of  units  allocated  to  IRL  (26)  and  DRL  (9).   The  minimum  total  annual 
equivalent  cost  is  $544,830.56. 

A  search  of  the  various  combinations  of  IRL  and  DRL  in  a  throw- 
away  age  is  executed  and  the  optimum  for  a  population  size  is  printed 
out.   In  each  population  the  number  of  IRL  channels  searched  will 
normally  be  between  1  and  12,  while  the  number  of  DRL  channels  can 
vary  from  0  through  6.   There  must  be  at  least  one  IRL  channel;  main- 
tenance must  be  performed  at  that  level  but  there  need  not  be  a  depot 
level  of  maintenance. 

Once  a  printout  is  available  the  minimum  cost  in  each  throw-away 
age  can  be  determined.   The  minimum  of  these  minima  is  the  global 
minimum.   In  the  printout  of  the  standard  example,  the  throw-away  age 
containing  the  minimum  is  10.   The  total  population  size  is  39,  and 
the  TAEC  is  $519,324.81. 

Figure  2  shows  the  minimum  TAEC  versus  throw-away  age,  this  plot 
shows  which  throw-away  age  has  the  minimum  cost.   The  associated  print- 
out must  be  used  to  determine  the  system  configuration  which  produces 
this  minimum  cost. 
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Figures  3  and  4  are  graphic  presentations  of  total  population  size 
versus  the  minimum  total  annual  equivalent  cost.   Each  throw-away  age 
is  represented  and  the  global  minimum  can  be  seen  in  Figure  4.   These 
figures  are  developed  from  the  standard  run  whose  output  is  included  in 
Appendix  B.3.   Two  figures  are  used  to  facilitate  reading  with  two 
different  vertical  scales  to  display  the  results. 

Discussion 

Analysis  of  the  output  in  Appendix  B.3  indicates  that  there  is 
considerable  information  available  to  the  designer  beyond  that  already 
discussed.   Should  he  have  the  option  to  use  a  range  above  the  minimum 
TAEC  for  the  system  to  be  procured,  there  are  several  system  con- 
figurations which  now  become  acceptable.   In  this  problem  if  the  limit 
was  plus  2  percent  of  the  minimum  cost,  he  would  then  look  for  all  com- 
binations with  a  cost  less  than  $530,400.   Scanning  the  output  data  it 
is  seen  that  there  are  14  options  within  the  2  percent  range  limit. 

Another  area  the  designer  may  wish  to  investigate  is  when  there 
is  a  maximum  amount  of  financial  resources  available  or  a  limit  on  the 
acceptable  annual  equivalent  cost.  He  can  quickly  ascertain  from  the 
output  if  there  is  an  acceptable  system  configuration.  He  may  need  to 
rerun  the  program  varying  the  system  parameters  until  he  attains  an 
acceptable  solution  or  decides  that  he  cannot  meet  the  demand  with  an 
acceptable  system  configuration. 
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CHAPTER  4 

ALGORITHM  SENSITIVITY  ANALYSIS 

This  chapter  presents  an  analysis  of  the  mathematical  model  and 
the  optimization  algorithm  developed  in  Chapter  2.   Several  system 
parameters  will  be  varied  and  their  effects  on  cost  will  be  investigated, 
This  analysis  will  allow  the  identification  of  critical  parameters. 

The  parameters  to  be  tested  will  have  their  values  specified  by 
the  designer.   These  parameters  will  be  assumed  to  be  logical  estimates 
and  will  be  substituted  into  the  standard  input  set.   The  original 
system  design  example  output  will  be  used  as  the  basis  for  all  compari- 
sons. The  variable  values  used  in  making  the  test  runs  will  be  listed 
with  each  test  discussion.   The  lines  in  the  standard  input  set, 
Appendix  B.2,  will  be  numbered  from  top  to  bottom,  one  through  15,  for 
purposes  of  reference  in  this  chapter. 

Reliability  and  Maintainability  Variation 

The  item  to  be  investigated  first  is  the  determination  of  the 
system  sensitivity  to  changes  in  MTBF  and  MTTR.  These  changed  factors 
result  in  either  the  increasing  or  decreasing  of  unit  availability. 
Availability  is  used  to  realign  the  unit  cost.   The  resultant  changes 
to  be  incorporated  in  the  new  input  set  are  given  below. 

It  is  assumed  that  for  each  one  percent  increase  in  availability 
the  initial  unit  cost  will  increase  by  $3,000.   The  cost  will  be  de- 
creased for  a  drop  in  availability  at  the  rate  of  $500  per  one  percent 
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change.   The  following  is  a  listing  of  the  system  input  data  used  in  the 
test  runs. 

First  Test  Run 

Line  4  MTTR:   same  as  standard  run 

Line  5  MTBF  (days):   9.1;  10.5;  12.0;  12.0;  12.0;  11.5; 

11.0;  10.5;  10.0;  9.5;  9.0;  8.5;  8.0;  7.0;  6.0. 
Line  6  Unit  cost  input:   $40,000;  $4,000  salvage  value;  15  year 

life;  Interest  rate  10%. 

Second  Test  Run 

Line  4  MTTR:   same  as  standard  run 

Line  5  MTBF  (days):   4.0;  5.0;  6.0;  6.0;  6.0;  6.5;  6.25; 

6.0;  5.5;  5.0;  4.5;  4.0;  3.5;  3.0;  2.0. 
Line  6  Unit  cost  input:   $23,100;  $2,310  salvage  value;  15  year 

life;  Interest  rate  10%. 

Third  Test  Run 

Line  4  MTTR  (days):   1.7;  1.5;  1.3;  1.25;  1.25;  1.25;  1.35; 

1.4;  1.45;  1.5;  1.6;  1.7;  1.8;  1.9;  2.0. 
Line  5  MTBF:   same  as  second  test  run 
Line  6  Unit  cost  input:   $21,800;  $2,180  salvage  value;  15  year 

life;  Interest  rate  10%. 
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Fourth  Test  Run 

Line  4  MTTR  (days):   1.3;  1.2;  1.1;  1.0;  0.9;  1.0;  1.1; 

1.2;  1.25;  1.3;  1.35;  1.4;  1.5;  1.6;  1.8. 
Line  5  MTBF  (days):   5.5;  6.9;  8.6;  8.6;  8.6;  8.0;  7.5; 

7.0;  7.0;  6.5;  6.5;  6.0;  6.0;  5.5;  4.5. 
Line  6  Unit  cost  input:   $29,500;  $2,950  salvage  value;  15  year 

life;  Interest  rate  10%. 

As  indicated  in  Figure  5,  it  is  more  economical  to  procure  and 
deploy  a  system  with  the  standard  system  design  parameters  than  any  of 
those  tested.   It  is  more  cost  effective  to  increase  the  unit  avail- 
ability, i.e.,  the  ratio  of  MTBF  and  MTBF  plus  MTTR,  than  to  decrease 
it.  A  negative  increment  of  the  availability  resulted  in  a  higher 
TAEC  than  the  same  increment  in  the  positive  direction  from  the 
standard  run  availability.   In  Figure  6  this  fact  is  also  shown. 

Percent  Population  Assigned  to  IRL  and  DRL 

In  this  test  the  percent  of  the  total  population  of  units  in  the 
system  allocated  to  IRL  and  DRL  is  varied.   The  interval  is  65  percent 
to  85  percent  inclusive,  in  five  percent  steps.  This  allows  a  repre- 
sentative cross-section  of  allocation  estimates  the  designer  might  use, 

First  Test  Run 

Line  2   IRL  percent  =  65 
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Second  Test  Run 

Line  2  IRL  percent  =  70 

Third  Test  Run 

Line  2  IRL  percent  =  80 

Fourth  Test  Run 

Line  2  IRL  percent  =  85 

As  depicted  in  Figure  7,  the  plot  of  the  TAEC  versus  throw-away 
age  indicates  that  there  is  very  little  variation  in  cost  when  compared 
to  the  standard  percentage,  75.   This  analysis  shows  that  the  population 
allocation  to  IRL  and  DRL  could  vary  between  65  percent  and  85  percent 
with  only  minor  changes  in  the  minimum  TAEC. 

Depreciation  Method 

The  next  parameter  modification  is  the  method  of  depreciation 
utilized  in  determining  the  book  value  of  the  units.   The  standard 
design  example  run  was  made  using  the  straight  line  method.   In  this 
analysis  the  other  three  methods  available  in  the  model  will  be  examined, 
First,  the  sum  of  years  digits;  second,  the  sinking  fund  method;  and 
thirdly  the  book  value  equal  to  zero  are  investigated. 

First  Test  Run 

Line  1  Depreciation  code  =  3  » 
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Second  Test  Run 

Line  1  Depreciation  code  =  4 

Third  Test  Run 

Line  1  Depreciation  code  =  5 

In  Figure  8  it  is  readily  apparent  that  the  two  modes,  book  value 
equal  zero  and  the  sum  of  years  digits  method  are  less  costly  during 
the  initial  ownership  years.  However,  they  both  approach  the  straight 
line  method  when  the  throw-away  age  approaches  the  maximum  life.   The 
sinking  fund  method  is  very  much  more  expensive  in  the  initial  owner- 
ship years.   It  too  approaches  the  other  methods  as  the  units  become 
older.   They  all  are  equal  at  the  end  of  the  life  cycle  of  15  years 
as  they  are  expected  to  be  in  this  model. 

Interest  Rate  Variation 

The  TAEC  will  be  varied  because  of  incremental  changes  in  the 
interest  rate.   The  effects  of  these  changes  will  be  investigated  in 
this  section. 

First  Test  Run 

Lines  6  through  12  Interest  rate  =  8% 

Second  Test  Run 

Lines  6  through  12  Interest  rate  =  9% 
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Third  Test  Run 

Lines  6  through  12  Interest  rate  =  11% 

Figure  9  shows  the  result  of  the  new  input  parameters  and  the 
computations  for  the  TAEC.   The  variation  is  equally  spaced  from  the 
standard  run  cost  for  a  particular  throw-away  age.   This  should  be 
expected.   The  difference  at  the  ten  year  minimum  cost  level  for  the 
standard  run  and  the  two  adjacent  percentages  is  approximately  $12,000 
above  and  below  the  standard  run  TAEC. 

Depot  Level  Repair  Time  Variation 

In  this,  the  final  sensitivity  analysis  test  to  be  run,  the  depot 
repair  time  is  varied.   This  test  will  identify  if  there  is  a  limit 
to  the  time  out-of -service  due  to  a  major  failure  which  results  in  the 
unit  being  sent  to  the  depot  facility  for  repair,  and  its  effect  on  the 
TAEC. 

First  Test  Run 

Line  3  DMTTR  (days):   5.5 

Second  Test  Run 

Line  3  DMTTR  (days):   12.5 

The  results  presented  in  Figure  10  shows  that  there  is  an  apprecia- 
ble increase  in  system  TAEC  when  the  mean-time-to  repair  is  increased 
by  800  percent.  However,  slight  changes  in  DMTTR,  up  to  about  5.5  days, 
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cause  only  moderate  increases  in  TAEC,  this  is  not  considered  an  ex- 
cessive variation.   This  fact  is  of  importance  to  the  designer  as  it 
allows  for  some  latitude  in  the  depot  facility  design  and  location  as 
well  as  the  amount  of  test  equipment  and  level  of  spares  that  must  be 
on-hand  for  effective  and  efficient  unit  repairs. 

Summary 

A  summary  of  the  results  of  the  sensitivity  tests  is  given  in 
Table  I.   The  table  lists  the  percent  change  in  the  parameter  varied, 
with  respect  to  the  standard  run  parameter,  and  the  resultant  change 
in  the  TAEC,  compared  with  the  minimum  TAEC. 
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TABLE  I.  Percent  Variability  of  TAEC  in  Sensitivity  Analysis 


Throw-away  Age  =10 

% (STANDARD) 

%(TAEC) 

I.  Reliability  and  Maintainability 

Run  No.  1 
Run  No.  2 
Run  No.  3 
Run  No.  4 

Std.  =  85.9% 

$519324.81 

+5.7 
-3.5 
-6.6 
+1.5 

+2.5 
+3.3 
+8.4 
+1.5 

II.  Percent  Population  IRL/DRL 

Run  No.  1 
Run  No.  2 
Run  No.  3 
Run  No.  4 

Std.  +  75% 

$519324.81 

-13.33 
-6.67 
+6.67 

+13.33 

+2.59 
+2.09 
-0.55 
-0.79 

III  Method  of  Depreciation 

Run  No.  1 
Run  No.  2 
Run  No.  3 

Std.  =  2 

$519324.81 



-4.76 
+3.01 
-2.66 

IV.  Interest  Rate  Variable 

Run  No.  1 
Run  No.  2 
Run  No.  3 

Std.  =  10% 

$519324.81 

-0.20 
-0.10 
+0.10 

-4.14 
-2.09 
+2.10 

V.  DMTTR  Variable 

Run  No.  1 
Run  No.  2 

Std.  +  2.5 

$519324.81 

+220 
+800 

+2.43 
+5.12 

CHAPTER  5 

ALGORITHM  CONSTRAINT  ANALYSIS 

This  chapter  will  extend  the  analysis  of  the  model  to  investigate 
the  application  of  constraints  on  the  design  variables  to  determine 
their  effects  on  system  effectiveness.   The  two  variables  that  will  be 
considered  are  the  throw-away  age  and  the  total  population  size.   The 
format  for  the  test  run  discussions  will  be  the  same  as  that  used  in 
Chapter  4. 

Throw-away  Age  Constrained 

The  throw-away  age  can  be  restricted  by  the  designer  or  by  manage- 
ment to  only  have  units  in  the  system  newer  than  a  prescribed  age.   In 
this  case  he  would  want  to  evaluate  the  effect  of  such  a  decision  on 
the  TAEC.   This  investigation  requires  four  of  the  data  input  lines 
be  modified.   They  are  as  listed  below. 

First  Test  Run 

Line  1  NNAGE  =  4 

Line  4  MTTR  (days):   1.5;  1.0;  1.25;  1.9. 

Line  5  MTBF  (days):   5.0;  8.0;  6.5;  4.0. 

Line  6  modify  life  of  units  to  4  years. 

Second  Test  Run 

Line  1  NNAGE  =  8 
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Line  4  MTTR  (days):  1.5;  1.1;  0.9;  1.2;  1.3;  1.4;  1.6;  1.9. 
Line  5  MTBF  (days):  5.0;  8.0;  8.0;  7.0;  6.5;  6.0;  5.5;  4.0. 
Line  6  modify  life  of  units  to  8  years 

Third  Test  Run 

Line  1  NNAGE  =11 

Line  4  MTTR  (days):   1.5;  1.3;  1.0;  0.9;  1.2; 

1.3;  1.4;  1.5;  1.6;  1.75;  1.9. 
Line  5  MTBF  (days):   5.0;  6.5;  8.0;  8.0;  7.0; 

6.5;  6.0;  5.5;  5.5;  5.0;  4.0. 
Line  6  modify  life  of  units  to  11  years 

In  Figure  11  the  program  output  is  plotted  as  previously;   TAEC 
versus  Throw-away  Age.   The  results  are  very  definite,  the  cost  incurred 
by  limiting  the  unit  life  is  very  high.  The  TAEC  does  follow  the  ex- 
pected trend  as  it  will  come  very  near  the  standard  design  run  as  it 
approaches  15  years.  At  eleven  years,  with  a  throw-away  age  of  10,  the 
difference  in  TAEC  is  plus  $50,000. 

Total  Population  Constrained 

There  may  also  be  some  physical  constraint  placed  on  the  scope  of 
the  system.   The  space  available  in  which  to  store  the  units  could  be 
a  very  serious  problem.   In  this  case  the  designer  would  have  to  limit 
the  number  of  units  to  be  procured  and  ascertain  what  effect  this  has 
on  the  TAEC.   In  this  problem  two  lines  of  input  data  need  to  be  modified. 
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First  Test  Run 

Line  1  XN  =  35 
Line  2  ENDXN  =  26 

Second  Test  Run 

Line  1  XN  =  30 
Line  2   ENDXN  =  21 

Figure  12  shows   the  expected  effect  of  a  constraint  and  also  an  un- 
expected effect.   In  the  case  of  the  second  test  run,  an  increase  in  the 
system  cost  is  exhibited  due  to  the  population  size  constraint.   However, 
the  first  test  run  shows  the  effect  of  rounding  the  factorial  and  com- 
binatorial calculations.   An  exhaustive  analysis  of  the  terms  in  the 
equations  described  in  Chapter  2,  and  contained  in  subroutine  UNA VAIL, 
lead  to  consideration  of  the  terms  involving  factorials  and  combinator- 
ials.   These  terms  are  rounded  and  as  a  result  the  accuracy  of  the 
equations  is  affected.   These  calculations  in  turn  result  in  the  varia- 
tion of  the  number  of  units  available.   This  number  is  relatively  small 
but  has  a  considerable  effect  on  the  cost  of  not  meeting  demand.   As 
an  example,  a  change  of  one-tenth  in  the  number  of  units  available  per 
day  below  the  established  minimum  required  in  the  algorithm  will  re- 
sult in  a  cost  of  $7,300  per  year.   In  the  case  analyzed,  the  unit 
availability  changed  sufficiently  to  incur  a  cost  on  the  system  in 
the  unconstrained  case  but  not  in  the  constrained  thereby  giving  the 
false  indication  of  a  lower  system  cost.  ' 
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Summary 

To  summarize,  the  constraint  analysis  shows  that  the  effect  of  con- 
straints on  system  cost  is  more  pronounced  and  they  have  greater  influence 
on  the  TAEC  than  system  parameter  variations. 


CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

As  a  result  of  this  study  it  is  concluded  that  a  multivariant 
system  for  a  finite  population  with  multiple  queue  channels,  a  two 
level  maintenance  capability,  and  optimization  of  the  system  to  meet 
a  demand  at  minimum  total  system  cost  can  be  achieved.   In  Chapter  3 
it  was  shown  that  utilization  of  the  designers  parameters  does  lead 
to  the  determination  of  a  global  minimum.  The  optimum  throw-away  age, 
the  number  of  intermediate  and  depot  level  facilities,  the  optimum  total 
population  size,  and  the  minimum  total  annual  equivalent  cost  for  that 
particular  system  configuration  can  be  specifically  determined  coinci- 
den tally. 

From  the  results  of  the  sensitivity  and  constraint  analyses  in 
Chapters  4  and  5  respectively,  it  can  be  concluded  that: 

1)  The  MTBF  and  MTTR  values  utilized  by  the  system  designer 
has  considerable  effect  on  the  overall  TAEC. 

2)  The  optimal  inherent  availability  of  the  units  to  be 
procured  should  be  approximately  86  percent  to  achieve 
the  most  effective  system  for  the  case  studied. 

3)  The  percent  allocation  of  the  failed  units  to  the 
intermediate  or  depot  level  of  repair  has  little 
effect  on  the  total  system  cost  in  the  range  analyzed. 
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4)  The  depot  MTTR  has  little  effect  on  the  TAEC  when  the 
increase  in  MTTR  is  less  than  200  percent.  Large  changes 
in  MTTR  do  increase  the  cost  by  more  than  five  percent. 

5)  The  interest  rate  and  the  depreciation  mode  used  in  the 
model  have  predictable  results.  All  the  results  were 
consistant  with  the  formulas  incorporated  in  the  model. 

6)  Limiting  the  throw-away  age  of  the  units  will  raise  the 
TAEC  of  the  system  if  the  limit  is  low  in  comparison  to 
the  unit's  original  design  life. 

7)  The  constraining  of  the  total  population  size,  combined 
with  the  system  parameters  estimated  and  utilized  in  this 
model  was  such  that  the  effects  of  computer  word  size  and 
rounding  were  dramatically  exhibited.   This  is  the  result 
of  fortunate  selection  of  system  parameters  and  variables. 
The  second  test  run  exhibits  the  normally  anticipated 
results. 

In  this  study,  the  concept  of  designing  the  total  system  from  the 
beginning  is  considered  to  be  the  critical  factor  in  achieving  the 
optimal  system  design.   The  demand  is  the  starting  point  in  arriving 
at  the  specification  of  a  system  which  will  optimally  meet  the  demand, 
The  total  population,  the  types  and  number  of  repair  facilities,  the 
reliability,  maintainability,  costs  of  unmet  demands,  initial  cost  of 
units  and  facilities,  and  the  necessary  support  and  logistics  were 
all  addressed  in  this  study. 
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The  mode  of  transformation  of  all  the  requirements  is  the  model 
and  its  associated  algorithm.   This  is  but  a  small  part  of  the  planning 
and  estimating  required  to  bring  a  system  into  being  efficiently  and 
effectively.   The  designer  can  now  write  the  specification  and  allocate 
the  requirements  to  the  various  engineering  and  manufacturing  specialists 
for  their  use  in  the  detailed  engineering  and  production  required  to 
bring  the  system  into  being. 

Recommendations 

There  are  several  areas  addressed  in  this  study  which  should  be 
analyzed  further.   They  generally  pertain  to  the  basic  system  design. 
A  few  of  these  are  as  follows: 

1)  Further  refinement  of  the  model  is  needed  in  the  areas  of 
support  and  personnel  costs.   The  estimates  made  are 
percentages  representative  of  systems  in  use  today.  More 
exact  methods  should  be  sought  to  delineate  these  costs 
more  accurately. 

2)  An  analysis  of  the  model  and  necessary  modifications  to 
the  algorithm  should  be  made  to  enable  more  detailed 
investigation  of  variations  in  demand  and  constraints  on 
the  number  of  IRL  and  DRL  channels  and  their  effect 

on  the  system  total  cost. 

3)  A  study  should  be  made  to  identify  the  parameters  and 
variables  effected  when  the  total  system  cost  is  limited 
to  some  maximum  amount.   This  will  require  a  substantial 
change  in  the  search  routine. 
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4)  An  analysis  of  the  out-of-operation  units  should  be 
made  to  establish  the  distribution  of  the  units  down. 
This  will  facilitate  establishing  a  better  equation 
for  the  unmet  demand  cost. 

5)  Continued  efforts  should  be  made  to  determine  the 
critical  system  parameters  early  in  the  system  design. 
The  algorithm  and  model  can  be  further  refined,  by 
updating  equations  and  reducing  the  number  of  assumptions 
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The  following  mathematical  development  is  contained  in  Chapter  III 
of  a  manuscript  entitled,  "Analysis  of  Queuing  Systems"  by  G.  K.  Bennett, 
J.  W.  Schmidt,  and  J.  A.  White  of  the  Department  of  Industrial  Engineering 
and  Operations  Research  at  Virginia  Polytechnic  Institute  and  State 
University. 

This  development  utilizes  the  classification  technique  introduced 
by  Kendall  (10)  .  The  basic  queue  developed  for  use  in  this  thesis  is  a 
Finite  Queue,  (M/M/c) : (GD/N/N)  is  the  type  which  is  developed  in  this 
Appendix.   In  order  of  occurrence,  the  letters  used  above  are  interpreted 
as  follows: 

M  -  denotes  Poisson  arrivals  to  the  repair  facility 

M  -  denotes  Poisson  repair  times 

c  -  denotes  the  number  of  repair  channels 
GD  -  indicates  the  units  are  served  according  to  any  general 
queuing  discipline 

N  -  is  the  maximum  number  of  units  allowed  in  the  repair 
system  (in  repair  and  waiting) 
and  N  -  is  the  finite  population  size. 

The  probability  mass  function  for  the  number  of  units  in  the  system 
is  as  follows: 

X  =  ..rPT,-,  =  expected  number  of  failures  per  unit  time 
MTBF       r  r 

y  =       =  expected  number  of  repairs  per  unit  time 

p  =  —   =  traffic  intensity,  system  utilization  factor,  or 
system  load  factor. 
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®    (cp)n       P  n=0,l,...,c+l 


P       =      < 


V  c!        V 


n=c,c+l, . . . ,N 


n>  N 


Solving  for  P     we  have: 


c-1    »  N    xt     ,  c  n 

i    0  <^n  +  i   0 


-1  -1 


11=0 


n=c 


c! 


To  determine  the  expected  number  of  units  in  repair,  L,  and  the  expected 
time  a  unit  waits  in  the  queue,  W  ,  let  the  following  equation  be  used; 
Peck  and  Hazlewood  (14) : 


X  = 


l/y 


1/y  +  1/X 


F  = 


1/X  +  1/u 
1/X  +  l/y  +  w 


F  represents  the  overall  efficiency  of  the  repair  system.   Then  based  on 
these  values  for  F  and  X,  the  following  equations  are  developed: 

L  =  N   [1  -  F  (1  -  X)] 


Wq  =   (1  -  F)/F  X 
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To  determine  W  in  terms  of  P  an  additional  value,  D,  delay  probability; 
is  defined  by  Peck  and  Hazlewood,  and  must  also  be  investigated  for  use 
in  solving  the  queuing  problem.   The  delay  probability  is  the  probability 
of  an  arriving  unit  not  being  able  to  proceed  directly  into  service  upon 
arrival  at  a  facility.   To  obtain  this  quantity  it  must  be  noticed  that 
it  is  no  longer  true  that  an  arrival  is  equally  likely  to  occur  in  any 
interval  of  time  At.   The  probability  of  an  arrival  during  At  is  dependent 
of  the  number  of  units  in  the  system.   Notice  that: 

P(unit  arrives  for  service  during  At) 

N 
=  £  P(unit  arrives  for  service  during  At/n  in  the  system)  P 

n 

N 

=  Y   (N  -  n)   P  X   At 
L  n 

n 

Now,  let  delay  probability  be  defined  as  follows: 

D  =  P(unit  arrives  for  service  during  At  and  has  to  wait) 
*  P(unit  arrives  for  service  during  At) 

N 

J   (N  -  n)  P   X  At 
**  n 

n=c 


thus,  D  = 


N 

I      (N  -  n)  Pn  X   At 


n  n 

n=0 


We  then  can  obtain  the  probability  mass  function  for  the  number  of  units 
in  the  system  ahead  of  the  arriving  unit  or  customer  by  noticing  that 

P(arrival/n  in  the  system) P 


P  (n  in  the  system/arriving  unit)  =  p (arrival) 


n 
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or  P 


n 


(N  -  n)  P 

n 

N 

I        (N  -  n)   P 
n=0 


In  designing  this  system  the  service  is  assumed  to  be  on  a  first 
come  first  served  basis.   Therefore  the  density  function  for  the  waiting 
time  can  be  expressed  as 


f(V 


N 


n=0 


f(t  /n)  P 

q      n 


where:   f(t  /n)  =  6(t  -0)  for  n=0,l,2, . . . ,c-l  and  f(t  /n)  is  the  prob- 
ability distribution  of  the  sum  n-c+1  independent  and  identically  dis- 
tributed random  variables  with  the  parameters  cu  where  n=c,c+l, . . . ,N-1. 
Thus,  f (t  /n)  is  a  gamma  distribution  with  parameter  n-c+1  and  cu  if  n  > 
c  or 


c-1 

Lr,         n 


fCtq>»< 


n=0 


N-l   ,  N n-c+1 
n=c 


(n-c)l 


>     \"° 


-cut 

t  n"c  e    q  P  ,   t  >  0 
q  n     q 


Since  E(t  ) 

q 


N 


n=c 


N  (t  /n)  P 
q     n 


we  see  that  for  n>  c 


-n-c+l> 


N-l 
q     ntc        {   cu  > 
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W 


N-l     ,,    (N-n)  P 
y   n-c+1         n 


n=c 


cy 


N 


I        (N-n)  P 

n=0        r 


N-l 


n=c 


n-c+1  ,..  N   N! 
(N-n) 


~c  n 

c  ,  P 


cy 


(N-n) !    c! 


W  = 

q 


c-l 


I      (N-n) 


n! 


N 


n=0 


nl(N-n)l 


(cp)n  +  I      (N-n) 


NI 


c  n 
c P 


n=c 


(N-n) !   c ! 


L  =  N 


1  - 


(1/y  +  1/X) 
(1/y  +  1/X  +  Wq) 


1   - 


1/y 
(l/y  +  iA)j 


APPENDIX  B 
B.l  Fortran  Program  Listing  of  System  Design  Model 
B.2  Input  Data  for  System  Design  Problem 
B.3  Output  Results  for  System  Design  Problem 
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LIST  OF  KEY  VARIABLES 


NNAGE  =  maximum  life  expectancy  of  a  unit 
NAGE   =  current  age  of  a  unit  for  calculation 
D      =  demand 

XN     =  total  number  of  units 
XNC    =  number  of  units  in  IRL 
XND    =  number  of  units  in  DRL 

NFCOST  =  number  of  cost  areas  for  facility  calculations 
ALEOC  ■  first  cost  of  establishing  IRL 
ALEOD  ■  first  cost  if  establishing  DRL 

ENDXN  =  population  size  for  termination  of  calculations 
AGEINC  =  increment  of  NAGE 

PERCNT  »  allocation  of  failed  units  to  IRL  and  DRL 
SUPPC  =  percent  allocated  for  recurring  support 
TRNPC  =  percent  allocated  for  training  and  init.  support 
P      =  initial  unit  cost 
SV     =  salvage  value  of  unit 
NP     =  life  expectancy  of  unit 
XINT   =  interest  (discount)  rate 

PFAC   =  initial  cost  of  IRL  facility  (PFACD  for  DRL) 
FFAC   =  salvage  value  of  IRL  facility  (FFACD  for  DRL) 
LFS    =  life  expectancy  for  IRL  (LFSD  for  DRL) 
DYCOST  =  admin.,  labor,  etc.  for  IRL 
DDCOST  =  admin.,  labor,  etc.  for  DRL 
CC     =  number  of  repair  channels  in  IRL 
DC     =  number  of  repair  channels  in  DRL 
ZN     =  daily  cost  of  not  meeting  a  demand 
TUAEC  =  total  unit  annual  equivalent  cost 
FAECN  =  facility  annual  equivalent  cost,  IRL 
FAECDN  =  facility  annual  equivalent  cost,  DRL 
WCNT   =  yearly  cost  of  not  meeting  a  demand 
TTAENC  =  total  cost  of  units  allocated  to  IRL 
TTAEND  =  total  cost  of  units  allocated  to  DRL 
TTAEST  =  cost  of  units  for  total  system 
FFAEST  =  cost  of  facilities  for  total  system 
STRN   =  training  and  initial  spares  cost 
RSUPN  =  recurring  support  cost 
TTAECN  =  grand  total  annual  equivalent  cost 
(program  search  based  on  this  term) 
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SELECTED  SENSITIVITY  AND  CONSTRAINT  STUDIES  FOR  A  FINITE 
POPULATION  QUEUING  SYSTEM  DEPLOYED  TO  MEET  A  DEMAND 

by 
Carl  K.  Hinger 

(ABSTRACT) 

The  purpose  of  this  thesis  is  to  study  the  sensitivity  of  certain 
parameters  in  a  conceptual  system  containing  a  finite  population  of 
units  deployed  to  meet  a  demand.   Selected  constraints  on  the  design 
variables  are  also  treated.   The  measure  of  effectiveness  in  this  system 
design  is  the  total  annual  equivalent  cost.   This  cost  is  the  summation 
of  the  total  population  unit  costs,  depot  and  intermediate  repair  fac- 
ility costs,  training  and  initial  support  costs,  recurring  support  costs, 
and  cost  of  not  meeting  a  demand  on  the  system. 

Three  separate  areas  are  addressed  in  this  study.   The  establishment 
of  a  standard  system  design,  the  study  of  the  variation  in  the  cost  as 
compared  to  the  standard  design,  and  the  effects  of  constraints  on  the 
system  cost  compared  to  the  cost  of  the  standard  design. 

A  special  computer  algorithm  was  developed  to  perform  a  multiple 
search  through  the  design  variables.   This  search  yields  the  minimum 
total  system  cost  and  a  resultant  system  configuration  to  be  specified 
by  the  designer  which  is  optimum.   The  algorithm  was  then  used  to  analyze 
the  effects  on  the  minimum  total  system  cost  of  variation  in  system 
parameters  and  constraints  on  design  variables. 
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